FINITE PROPAGATION SPEED FOR FIRST ORDER SYSTEMS 
AND HUYGENS' PRINCIPLE FOR HYPERBOLIC EQUATIONS 
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Abstract. We prove that strongly continuous groups generated by first order 
systems on Riemannian manifolds have finite propagation speed. Our procedure 
provides a new direct proof for self-adjoint systems, and allows an extension to 
operators on metric measure spaces. As an application, we present a new approach 
■ to the weak Huygens' principle for second order hyperbolic equations. 

a 

a. 

• 1. Introduction 

For a self-adjoint first order differential operator D acting on a space L 2 (V). 



Oh 

< 



where V is a vector bundle over a complete Riemannian manifold M , it is known 
that the unitary operators e ltD act with finite propagation speed. Indeed, if the 
principal symbol satisfies 

(1.1) Wd{x,0\ < k\£\ Wx G M, V£ G TIM 

for some positive number k, and if sppt(w) C K C M , then 
(S| : sppt(e itD u) C K K \ t \ := {x G M; dist(x,iv") < re|t|} Vi G M, 

>: 

that is, e has finite propagation speed kd < k. See, for instance, Proposi- 
tions 10.2.11 and 10.3.1 in [T3]. 

Noting that o~d{x, Vt)(x))u(x) = [D,rjl]u(x) for all bounded real-valued C 1 func- 
tions rj on M, we see that condition (II. ip implies that the commutator [rjI,D] is a 
O ! multiplication operator which satisfies the bound 



OO 



X 



(1.2) \\[r]I , D]u\\ 2 < kWVtjWooWu]^ Vn G Dom(£>) C Dom(Drj), 

where Dom(D) denotes the domain of D and [r)I,D]u = nDu — D(rju). 

On stating the result in terms of commutators, we can remove the differentiability 
assumptions altogether and consider operators D defined on metric measure spaces 
instead. Our aim is to also weaken the self-adjointness condition on D to the 
requirement that iD generates a Cq group (e ltD ) t m with 

(1.3) \\e UD u\\ 2 < ce w|t| ||u|| 2 Vt G R, Vw G L 2 (V) 

for some c > 1 and > 0. (When D is self-adjoint, condition (II. 3p holds with 
c = 1 and w = 0.) At the same time we could replace the L 2 space by an LP space 
if we wished. This requires a new proof of finite propagation speed. Let us state the 
result here for Riemannian manifolds. 
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Theorem 1.1. Let D be a first order differential operator which acts on a space 
L 2 (V) , where V is a complex vector bundle with a Hermitian metric, over a separable 
Riemannian manifold M . Suppose that iD generates a C group (e ltD ) t& m. satisfying 
( 11.3ft and that the commutators of D with bounded real-valued C°° functions rj on M 
satisfy ( II. 2ft . Then the group [e ltD )t£R has finite propagation speed kd<ck. 

In particular, the operator D acting on L 2 (V) could denote a first order system 
acting on L 2 (R™, C ) for some positive integers n and N , or on L 2 (Q, C ) where 
Q is an open subset of M. n . We remark that the constants in (11.21) and (II. 3p could 
be with respect to another norm on L 2 (V) equivalent to the standard one. We shall 
return to this point. 

In Section [3J we prove our main result, Theorem 13.11 which is a generalisation 
of Theorem 11.11 to metric measure spaces. The proof utilises a higher-commutator 
technique introduced by Mcintosh and Nahmod in Section 2 of [IE], and used to 
derive off-diagonal estimates, otherwise known as Davies-Gaffney estimates, by Ax- 
elsson, Keith and Mcintosh in Proposition 5.2 of [5], and by Carbonaro, Mcintosh 
and Morris in Lemma 5.3 of [S]. The proof also simplifies the argument based on 
energy estimates that is known for self-adjoint operators. 

A weak Huygens' principle for second order hyperbolic equations on IR n is proved 
as an application in Section [31 Homogeneous and inhomogeneous hyperbolic equa- 
tions are treated separately. The homogeneous version in Theorem 15.21 only requires 
the finite propagation speed result for self-adjoint first order systems. The inhomo- 
geneous version in Theorem 15.61 however, requires the generality of Theorem 13.11 
These results are achieved by introducing a first order elliptic system BD, where 
D is a first order constant coefficient system, and B is a multiplication operator, 
such that the second order hyperbolic equation contains a component of the system 
(BD) 2 . This approach is motivated by the work of Auscher, Mcintosh and Nah- 
mod @], and of Axelsson, Keith and Mcintosh [Sj, in which the solution of the 
Kato square root problem for second order elliptic operators is reduced to proving 
quadratic estimates for related first order elliptic systems. See also the survey by 
Auscher, Axelsson and Mcintosh [3]. 

2. Notation 

Henceforth, M denotes a metric measure space with a metric d(x, y) and a ct -finite 
Borel measure \x . In particular, M could be a Riemannian manifold as in the Intro- 
duction, or an open subset of M n with Euclidean distance and Lebesgue measure. 
If K,K C M and x G M, then d(x,K) := mf{d(x,y) ; y G K} and d(K, K) := 
M{d(x,y) ; x G K,y G K} . We define, for r > 0, K T := {x e M ; d(x,K) < r}. 
By Lip(M) we mean the space of all bounded real-valued functions rj on M with 
finite Lipschitz norm 

iuii = sup 

IK/ II Lip hup d , s 

Whenever K C M and a > 0, the real- valued function rjK,a defined by 

(2.1) f]K,a{ x ) '■= max{l — ad(x, K) , 0} 

belongs to Lip(M) with 1 1 »7jc,q: 1 1 Lip < a , an d sppt(r)K, a ) C Ky a . 

A vector bundle V over M refers to a complex vector bundle 7r : V — > M equipped 
with a Hermitian metric (•, -) x that depends continuously on x G M . The examples 
to be presented in Section [4] will be trivial bundles M x C m with inner product 
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(CjOx — SjO^j- F° r every vector bundle V, there are naturally denned Banach 
spaces L P (V), 1 < p < oo, of measurable sections. In particular, L 2 (V) denotes the 
Hilbert space of square integrable sections of V with the inner product (u , v ) := 
J M (u(x) , v (x)) x d/i(x) . In the case of the trivial bundle M x C m , these are denoted 
as usual by L p (M,C m )- 

The Banach algebra of all bounded linear operators on a Banach space X is 
denoted by C(X). Given A G L°°(M, C(C m )) , the same symbol A is also used to 
denote the multiplication operator on L p (M,C m ) defined by u \— > Au. Note that 
H^MIp < HAIIoollullp. Multiplication operators on L P (V) are defined in the natural 
way. For any function rj G Lip(M), the multiplication operator rjl : L P (V) — > L P (V) 
is defined by (rjl)u(x) := rj(x)u(x) for all u G L P (V) and /i-almost all x G M . This 
is a multiplication operator by virtue of the facts that 77 is bounded and continuous, 
and fi is a Borel measure. The commutator [A, T] of a multiplication operator A 
with a (possibly unbounded) operator T in L P (V) with domain Dom(T) is defined 
by [A,T]u = ATu - TAu provided u G Dom(T) n Dom(7M). 

Given an operator D in L P (V) (1 < p < 00), we say that i.D generates a Cq 
group (V(t))t £ R provided t H- V(£) is a strongly continuous mapping from R to 
£(LP(V)) with V(s + t) = V(s)V(t), 7(0) = J (the identity map on L"(V)) and 
|7(t)w| t=0 = iDu for all u G Dom(D) = {m G LP(V) ; |7(t)w| t=0 exists in L P {V)}. 
We write V(t) = e itD . Such a group automatically has dense domain Dom(Z^), 
and satisfies an estimate of the form ||e liD || < ce w '*' for some c > 1 and > 0. 
An introduction to the theory of strongly continuous groups can be found in, for 
instance, [H] or [TT]. We remark that, when D is self-adjoint in L 2 (V), Stone's 
Theorem guarantees that the operators e ttD are unitary, so iD generates a Co 
group with c = 1 and lo — . 

The group (e JtD )t G R is said to have finite propagation speed when there exists a 
finite constant k > 0, such that for all u G L P (V) satisfying sppt(n) C K C M , and 
all £ G R, it holds that sppt(e rfI? M) C K K \ t \. The propagation speed k^, is defined to 
be the least such k. 

3. The Main Result 

The following theorem is the main result of the paper. Theorem 11.11 is proved as 
a special case at the end of the section. 

Theorem 3.1. Let V denote a complex vector bundle over a metric measure space M 
and let 1 < p < 00. Suppose that D : Dom(D) C L P (V) — > L P (V) is a linear opera- 
tor with the following properties: 

(1) there exist finite constants c > 1 and u > such that iD generates a Cq 
group (e itD )tm in L P (V) with He*"^ < ce w l*l||n|| p Vt G R,Vu G L P (V); 

(2) there exists a finite constant k > such that for all 77 G Lip(M), one /ias 
77U G Dom(D) and || [n/, -D]u|| p < «||?7||Lip||w||p and [rjl, [rjI,D]]u = for all 
u G Dom(£)) . 

T/ien £/ie ^rowp (e** i: ')t e ]R /ias /im£e propagation speed k d < ck. 

Remark 3.2. As the commutator [r]I,D} is bounded on the dense domain Dom(D) , 
it extends uniquely to an operator (denoted by the same symbol) [rjl, D] G C(L P (V)) 
with the same bound. 
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Remark 3.3. The theorem remains true when the norm on L P (V) is replaced by 
another equivalent norm. This is used later in Case II of Section HI 

Remark 3.4. The results remain true in Bochner spaces L 2 (V) when the fibres of V 
are infinite dimensional Hilbert spaces. 

For completeness, we prove a known formula for the commutator [ql, e ttD ] . 

Lemma 3.5. Under the hypotheses of Theorem \3.1\ the following holds: 

[riI,e itD ]u = it [ e istD [r]I,D]e i{1 - s)tD uds Wt G R, Wu G L P {V). 
Jo 

Proof. It suffices to verify the expression when u G Dom(D). The property that 
{e ltD )t£R is a Co group then guarantees that e ltD u G Dom(Z^) with derivative 
j- t (e ltD u) = iDe ltD u = ie ltD Du for all t G M. The property that rju G Dom(D) 
then implies that e lstD ^j'j e l ( 1 - s ) tD u [ s different iable with respect to s, with 

^ stD {r,I)e^ tD u) = (e^faiy' 1 -^)'^) + e istD {[r,I)e^ l ~ ^ tD u)\s) 

= -ite istD [r]I,D]e i{1 - s)tD u 

for all s G R. This version of the chain rule can be found in, for instance, 
Lemma B.16 in [TT]. Using the fundamental theorem of calculus, we then have 

[rjl, e itD ]u = - I (e^iriiy^-^u) ds = it f e istD [ V I, D]e i{1 - s)tD u ds 
Jo ds \ / Jo 

as required. □ 

Proof of Theorem \3.1[ Given t G M and u G L P (V) with sppt(w) C K C M, our 
aim is to prove that sppt(e l * D w) C K CK \ t \. To do this, it suffices to prove that 
[e ltD u,v) = for all v G L p> (V) with d(sppt(u), K) > cn\t\ (where p' = ^j). Let 
us fix K,t,u,v, and choose a > such that c«|t| < 1/a < d(sppt(u), K) . On 
defining the cut-off function r) := r)K, a G Lip(M) as in (12. ip . we have rju — u, 
rjv — and cn\t\ IMIup < c«|t|a < 1. 

To simplify the computations, set 5 : C(L P (V)) — >■ £(L P (V)) to be the derivation 
defined by S(S) = [r)I,S] for all S G £(L P (V)), and adopt the convention that 
S°(S) := S. We see that 

(3.1) (e itD u,v) = (J tD ri n u,v) = -{8{j tD )r) {n - l) u,v) =...= (-l) n (6 n (e itD )u,v) 

for all n G N. The derivation formula 5(ST) = 5(S)T + S5(T) is readily verified for 
any S, T G £(L P (V)). Using Lemma [3.51 and the fact, given by property (j2J) of D, 
that S([t)I,D]) = [r]I, [i]I,D]] = 0, we then obtain 

„x m 

(3.2) 5 m+l {e ltD )u = it ^(^)5 m - fe (e iS ' D )[r//,D]5 fe (e l(1 - s) * D )M^ 

^° fc=0 

for all m e N := N U {0}, where the binomial coefficient (™) := ]^0zm ■ We now 
prove by induction that 

(3.3) \\S n (e UD )\\ < (c\t\ \\[r]I,D]\\) n ce"W 
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for all n E Nq. For n — , this is given by property (JT]) of D . Now let m G N and 
suppose that (13.31) holds for all integers n < m. We then use (13.21) to obtain 

\\5 m+1 (e itD )\\ < \t\ / ^(™)||5 m - fc (e is * D )|| || [77/, D] || \\5 k (e^- s)tD ) \\ ds 

Jo k=0 

-1 m 

< (c |t| || [rjI,D] || ) m + 1 ce "l'l / 

= (c|f|||[^,D]||)^ 1 ce u '* 1 / {s + {l-s)) m ds 

Jo 

= (c\t\\\[r ] I,D]\\r +1 ce^ . 

This proves ( 13. 3ft for all n G No. 

Therefore, using the estimate ( 13 .3p in ( 13. ip . together with property (EJ, we obtain 

\{e itD u,v)\ < (cre|f|||r7|| L ip) n ce w| * l ||w||p||'y|| p / < {cn\t\a) n ce u ^\\u\\ p \\v\\ pl 

for all n G N . We have c«|t|a < 1, so (e ltD u,v) = as required. □ 

Remark 3.6. In fact we have proved the stronger statement that 

s PP t(e UD u) C K t (D) C K CK]tl , 

where K t (D) = n{sppt(r?) ; 77 G Lip(M) , 77 = 1 on K , c\t\ || [rjl, D] < 1}. 

For example, if M = R> and ^ doe s not appear in D, then there is no propa- 
gation in the x\ direction. See the recent paper of Cowling and Martini [9] for some 
related results. 

We conclude the section by proving Theorem 11.11 

Proof of Theorem li.il From what was explained in the introduction, together with 
the fact that [77/, [rjl, D]] = when [77/, D] is a multiplication operator, Lemma l3~5l 
also holds under the hypotheses of Theorem 11.11 We then consider e > and 
repeat the proof of Theorem 13. 11 replacing the cut-off function r)K, a G Lip(M) with 
a [0,1] -valued function f] G C°°(M) satisfying f] = 1 on K, sppt(ry) C Ki/ a , and 
llV^Hoo < (1 + e) 1 1 rjK,a 1 1 Lip- These approximations exist because the Riemannian 
manifold M is separable (see, for instance, Corollary 3 in [7j). The result follows 
because e > can be chosen arbitrarily. □ 

Remark 3.7. The only reason that the Riemannian manifold M was required to be 
separable in Theorem 11.11 was so that we could construct smooth approximations to 
Lipschitz functions in the proof above. Indeed, Theorem 13.11 provides an analogous 
result without requiring separability. 

4. Some Special Cases 

A typical example of a first order system is the Hodge-Dirac operator D = 5 + 5* 
acting on L 2 (A(M)) when M is a complete Riemannian manifold. For this operator, 
the group (e ltD )t£R has finite propagation speed 1. We shall consider the case when 
the manifold is M" or an open subset thereof, and restrict attention to the leading 
components of the Hodge-Dirac operator (the components acting between scalar- 
valued functions and vector fields), and perturbations thereof. For this purpose, 
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when Q is an open subset of W 1 , we require the Sobolev space W 1,2 (Vt) consisting 
of all / in L 2 (Q) with generalised derivatives satisfying 



liyi.2(n) •— \\j \\L 2 {n) + IIWI 
where V/ = (djf) j=1) ... >n . 

Case I. Let M = M. n (n G N), and let D denote the self-adjoint operator 

W 1,2 (R n ) L 2 (W l ) L 2 (R n ) 



lL 2 (tt) < 00 : 



D 




V 



div 




C © -> 

L 2 (R n ,C n ) L 2 ( 



Dom(div) L 2 (M n ,C r 

where V : / \-> (djf)j has domain W 1 ' 2 ^" 1 ), and div = -V* : (uj)j H- ^jdjUj 
has domain {m G L 2 (M n ,C n ) ; divw G L 2 (M ra )}. 

It follows from known results that (e^)^ has finite propagation speed 1. It is 
also a consequence of Theorem 13.11 with c = 1 and u> = because Z) is self-adjoint, 
and with k — 1 because (using HV77H00 = | Ml Lip) we have 



[r,I, D] 



u 



3)3 



< d 3V)jf 



< 



WWup 



f 



3)3 



Vr] G Lip(M n ). 



Case II. Let D denote the operator in Case I, and consider the perturbed operator 
BD with the same domain, where B G L°°(W n , £(C 1+n )) satisfies (B(x)(, () > \\(\ 2 
for a.e. and all £ G C 1+n , for some A > 0. The multiplication operator B is a 

strictly positive self-adjoint operator in L 2 (R™, C 1+n ) , since (Bu, u) > A||w|| 2 . Hence 
||5 1//2 || = || B || l l 2 = || I? Hoc, 1 / 2 and B^ 1 , B~ 1 / 2 both exist as bounded operators. 
Using these facts, we find that BD is self-adjoint in L 2 (IR n , C 1+n ) under the 



inner product (u,v)b '■— (B u,v), whose associated norm \\u\\ B 
equivalent to ||u||. So iBD generates a C group (e ttBD ) te $i in L 2 (I 



W tBD u\\ < 



151 



1/2 11 itBD 



u\\b<\\B 



1 1/2 11 I, 

00 \\ U \\B 



< X-^WBl 



1/2 



llB-^uw is 
C 1+n ) with 

HI • 



The commutator [rjI,BD] = B[rjI,D] is a multiplication operator which satisfies 
\\[r]I, BD]u\\ < WBWooWr/WupWuW , so by Theorem I3.1[ the group (e ltBD ) t £R has finite 
propagation speed kbd < A~ 1 / 2 ||-E>|| 00 3//2 . 

Actually, this can be improved. As noted in Remark 13.31 the equivalent norm \\u\\b 
on L 2 (M. n , C 1+n ) can be used in the proof of Theorem 13. 11 The operator BD is self- 
adjoint in this norm, and \\[r]I, BD)u\\b = \\B^ l / 2 B[r]I , D]u\\ < ||-E>||oo|M|Li P |M|,B, 
so we conclude that (e ttBD ) t £M has finite propagation speed kbd < Halloo- 

Case III. Now we allow inhomogeneous terms, and consider an operator of the 
form BD acting on an open subset fl C IR n . Suppose that V is a closed subspace of 
W 1,2 (Q) which contains C£° (the C°° functions with compact support), and which 
has the property that T]f G V for all 77 G Lip(fi) and / G V. For example, the 
space W 1,2 (Q) itself has this property, as does W ' 2 (Q) (the closure of C£°(0) in 
W 1,2 (Q)). (This last statement follows from the facts that, given 77 G Lip(O), then 
rjf G Wo' 2 (fl) for all / G C c °°(fi) and rjl G £(W l > 2 (Q)) .) 

Define V v ■ V C L 2 (Q) ->■ L 2 {Q,C n ) by V v f = {d/f)j, and set divy = -Vy*. 
That is, divyu = divw for all u G Dom(divy) = {u G L 2 (fi,C™); divw G L 2 (f2) 
and (— V/,tt) = (/,divit) V/ G V"}. In particular, we have 771* G Dom(divy) for 
all 77 G Lip(O) and m G Dom(divy). 
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Define the self-adjoint operator 









V 


L 2 (n) 


L 2 (n) 







/ 


- divy 


© 


© 


© 


D = 


I 








: L 2 (n) 


C L 2 (9) 


-> L 2 {9) 




Vy 








© 


© 


© 








Dom(divy) 


L 2 (n,c n ) 


L 2 (n,c 



and the multiplication operator 



B 



a 







^00 (A ok ) 

O^o) (A jk ) 



e L°°(fi,£(C 2+n )) 



with a(x) > A > and Y%,k=i A 3k( x )(kCj > MC\ 2 for a.e. x G ft and all C G C n , 
for some A > 0. 

The operator BD satisfies (11.21) as in Case II. If A, and hence B , were positive self- 
adjoint, then iBD would generate a Co group as before, but we have only assumed 
this for the matrix- valued function (Ajk) with j, k G {1, . . . ,n} . We remedy this by 
writing 

a —a div 

BD A 00 + J2 k A ok d k =BD + C 
(A ]0 + E k A jk d k ) 



c 



where B 



ReA 00 + a (A ok ) 
(Aoj) (A jk ) 





i Im A 00 — a 
(^o-3o~) 



and a > 



is chosen large enough so that (B(x)£,Q > f|C| 2 fo r a - e - x E Q and all C G 
C 2+n . As in Case II, we see that BD is self-adjoint in L 2 (Q,C 2+n ) under the inner 
product (u , v) := (B~ 1 u , v ) , so i_BD generates a C group with ||e 



itBD 



u\\b < 



Lemma 14.11 below then allows us to deduce that iBD 
group (e ltBD ) tm with 



iBD + iC generates a C 



\e UBD u\ 



< \\B\ 



V2, 



Jt(BD+C) 



U \\B < 



\B\ 



< ce' 



^1*1 



\u\ 



for some finite c > 1 and u> > 0. By Theorem 13.11 we conclude that (e ltBD ) 
finite propagation speed. 



has 



Lemma 4.1. Let X be a Banach space, and suppose that T : Dom(T) C X — > X is 

a linear operator that generates a Cq group (e tT ) t £R in X satisfying \\e tT \\ < ce^'*' 
for some c > 1 and u > 0. If B £ t/ien t/ie sum T + £? on Dom(T) generates 

a C group (e t{T+B) ) t m on X satisfying ||e* (T+B) || < ce {uJ+c ^ Bm . 

Lemma 14.11 is a well known result based on the work of Phillips in [T7j . The proof 
for semigroups in Theorem III. 1.3 of [TT] can be extended to give the above result. 



5. Weak Huygens' Principle 

In this section, we apply Theorem 13.11 to prove a weak Hugyens' principle for 
second order hyperbolic equations. For motivation, we start with the wave equation 
on R™. A homogeneous equation with bounded measurable coefficients is treated 
next, followed by an inhomogeneous version on a domain Q C M. n . 
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In Case I, we considered the operator D 



-div 
V 



L 2 (R n ) 
in © , which 

L 2 {R n ,C n ) 

is self-adjoint, and noted that iD generates a Co group (e ltD )teR with finite propa- 
gation speed 1. Consequently, the cosine family cos(tD) = \{e %tD + e~ ltD ) (t > 0) 
also has finite propagation speed 1, where this is defined in the obvious way. Note 
that the cosine operators, being even functions of D, satisfy cos(tD) = cos(t\/D 2 ~), 

where D 2 = ^ y'div an d A = div V denotes the Laplacian operator with 

domain Dom(A) = {/ G W 1,2 (R n ) ; V/ G Dom(div)}. On restricting attention to 
the first component, we deduce that the cosine family (cos(t\/— A)) t > has finite 
propagation speed 1. This is at the heart of the weak Huygens' principle for the 
wave equation: 

Theorem 5.1. If f e W 1 ' 2 ^) , g G L 2 (R n ) with sppt(/) U sppt(#) C K C R n , 
then the solution 

F(t) = cos(tV^A)f + / cos(s\^A)gds 

Jo 

of the Cauchy problem 

§F(t) - AF(t) = (t > 0) 
KmF(t) = f 

has support sppt(F(t)) C K t . 

This result is very well known. The solution F belongs to C^R" 1 ", L 2 (R n )) n 
C°(R + ,W 1 ' 2 (R n )). There is a considerable literature on the wave equation, so we 
shall not proceed further with statements of uniqueness, energy estimates, etc. 

We turn now to the corresponding result for homogeneous equations with L°° 
coefficients. Let L = — adivAV = -aE^=i^)AA» where a G L°°(IR n ) and A = 
(A jk ) G L°°(R n ,£(C n )) with a(x) > A > and (A(x)(,() > MC? f o r a.e. x G R n 
and all ( e C n . Here L : Dom(L) C L 2 (R n ) -> L 2 {R n ) with Dom(L) = {/ G 
W 1 ' 2 ^) ; AVf G Dom(div)}. 

Theorem 5.2. If f E W 1 ' 2 ^), g G L 2 (M n ) sppt(/) U sppt(#) C K C R n , 
then the solution 

F(t) =cos(tv / L)/+ / cos(s\^L)g ds G C 1 (IR + , L 2 (R")) fl C°(1R + , W^ 1,2 (1R")) 
o/ t/ie Cauchy problem 

g>F(t) + LF(t) = (t > 0) 
limF(0) = / 

lim|F(0)=„ 

/ias support sppt(F(t)) C K at , where CY '. — ( || Q,|| oo || A\\ oq) ^' 2 . 
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Proof. Apply Case II with B 



- 

P 

(3A 



and j3 



l|A|| c 



1/2 



Then 



and (BDf 



BD 

L 




a 

/3 







' - div " 




-J div 


_ 


(3 A _ 




V 




I3AV 




L 



C group (e 
cos(tBD) = \ 



tBD 



with L as above and L = — AVadiv. 



) t£ K has finite propagation speed k^d < ||-B||oc 
ltBD ) i it is clear that the cosine family 



> e itBD + e 



From Case II, the 
On defining 



a . 



{cos(tBD)) t > = (cos(t^(BDy)) 



cos(tvT) 
cos(t 



t>0 



has the same bound on its propagation speed. It follows that the first component 
(cos(£\/L)) f >o, acting on L 2 (R™), has the same bound a on its propagation speed, 
as required. □ 

Remark 5.3. The operators cos(tBD) agree with those defined in the standard func- 
tional calculus of the self-adjoint operator BD acting on L 2 (R n , C 1+n ) with inner 
product (u, v)b = (B~ l u,v). Their properties are those of cosine functions as pre- 
sented, for example, in [U [2]. Note that the use of the square root sign is purely 
a symbolism to express the fact that the cosine operators cos(tBD) are diagonal. 
The function cos(t\/i) is an analytic function of z. 

Remark 5.4. When a = 1, then L is the self-adjoint operator in L 2 (IR n ) associated 
with the sesquilinear form J a ■ W 1 ' 2 (W l ) x iy 1,2 (]R n ) — >■ C defined by 

„ n 

JA(f,g)= / V] A jk (x)(d k f(x))d j g(x)dx 
Jw> j,k=i 

for all /, g G W 1,2 (M. n ). See, for example, Chapter IV of [Hj. The weak Huygens' 
principle is well known for these operators. See, for example, [IE]. Degenerate 
elliptic operators are also treated in [10] . 

Remark 5.5. Our methods work also for systems where the functions are C^-valued 
for some iV £ N, and A £ L°°(C nN ) satisfies Garding's inequality: 

ML f) > A||V/||i 2(RW)CnJV) V/ £ W 1,2 (R", C^) . 

The proof that iBD generates a C group needs to be modified as follows. 

The positivity condition on B £ L°°(M n , C(C N+nN )) is weakened to (BDu,Du) > 
\\\u\\ 2 for all u £ Dom(D). Then, following |3j, L 2 (R n ,C N+nN ) = N(£>) © R(BD) 
with corresponding projections Pn and Pr. (Here N(D) denotes the nullspace of 
D and R(BD) denotes the range of BD) The projections, which are typically non- 
orthogonal, commute with BD. Moreover B : R(D) — > R(BD) has a bounded 
inverse B^ 1 : R(BD) R(D) . The operator BD is self-adjoint in L 2 (R n , £ N+nN ) 
under the inner product {u ,v)b '■= (Pn^ , Pn^) + (B^Pru , Prv) , whose associated 
norm is equivalent to ||n||. Thus iBD generates a Co group. Proceed as in 

Case II, though note that the bound on kbd is not the same as before. 

Finally we consider inhomogeneous hyperbolic equations on a domain Q C M". 
As in Case III, suppose that V is a closed subspace of W /1,2 (fi) which contains 
C£°(Q), and which has the property that rjf £ V for all 77 £ Lip(f2) and / £ V. 
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with a(x) > A > and 



Suppose also that a,Aj k G L°°(IR n ) (j, k = 0, ...,n) 
Elk=i A 3kCkCj > A|C| 2 for a.e. xeHand all C G C". 

Define Vy : V C L 2 (ft) ->• L 2 (ft,C n ) by Vy/ = (djf)j, and set div v 



That is, divy u = divu for all u G Dom(divy) 
and (-V/,«) = (/,div«) WfeV}. 
Define the operator L in L 2 (Q) by 



{u e L 2 (nx n ); divu e L 2 (tt) 



Lu = —a djAjkdkU — a djAj u + a A kdkU + aA 00 u 



j,k=i 



3=1 



k=l 



[ I ~ diVy ] 



u 



A m (A ok ) I 
(A j0 ) (A jk ) \ [ Vy 

with Dom(L) = {u E V ; (J2k=i AjkdkU + Aj u)j G Dom(divy)}. 

Theorem 5.6. If f G V , g G L 2 (Vt) with sppt(/) U sppt(p) C K C Vt, then the 
solution 

(5.1) F(t) = cos(ty/L)f + / cos(sVL)g ds 

Jo 

of the Cauchy problem 



jgjF(t) + LF(t) = (t > 0) 



limF(0) = / 
hm|F(0)=^ 



has support sppt(F(t)) C K & t '■= {x G dist(x,K) < at} for some finite a. 

Proof. With B and D specified as in Case III, we have shown that iBD generates 
a Co group (e itBD ) t m with finite propagation speed that we now denote by a. The 
cosine family cos(tBD) = \{e itBD + e ' itBD ) = cos(t^(BD) 2 ) (t > 0) has the same 
propagation speed. Noting that 

(BD) 2 = 



L 

A)0 (Lok) 

o (L j0 ) (L jk ) 



for suitable operators Lj k , it follows that the first component (cos(ty/L)) t > , acting 
on L 2 (IR n ), has finite propagation speed bounded by a. The result follows. □ 

Remark 5.7. The use of the square root symbol is again purely symbolic, as cos(tV^) 
is really just the leading component of cos(tBD) . It is not the case in general that 
an operator a/L is defined. See Remark 15.31 In the language of cosine families, it 
is said that (BD) 2 , and hence L, are generators of their respective cosine families. 
In the usual treatment of cosine families associated with L, one adds a positive 
constant to A 00 if necessary to ensure that Re Ja(/, /) > || V/|| 2 + |M||> an d notes 
that the numerical range of A is contained in a parabola. This ensures that L 



generates a cosine family [12], and that ||vL/||2 ~ ||V/||; 



15] . See also pQ. 



Our treatment is consistent with this approach, but we do not need to apply it, 
as we use the more straightforward fact that, since iBD generates a Co group, 
then (BD) 2 generates a cosine family, and so then does the restriction to its first 
component L in L 2 (R"). 
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Remark 5.8. Given / G V and g G L 2 (R"), the solution F of the Cauchy problem 
satisfies F G L 2 (W 1 )) n C°(R+, V) with 



+ IIVF 



|^F(t)|| 2 <c(l + t)e' 



it 



+ l|v/|| 2 + N| 2 } 



for all t > and some constant c. Let us prove this using what we know about the 
operator BD in Case III. The constant Co is the same one as there. 

When u G L 2 (R n , C 2+n ) , then e itBD u G C°(R, L 2 (R™, C 2+n )) , and when u G 
Dom(D) = Dom(BD) , then e ifcBD u G C^R, L 2 (R™, C 2+n )) n C°(R, Dom(D)) . Hence 
a similar statement holds for cos(tBD). On restricting to the first component, we 
deduce that cos{ty/I)f G C\R + , L 2 {R n )) n C°(R + ,y) with 

|| cos(tVL)/|| 2 + ||Vcos(tVX)/|| 2 + ||4cos(*/L)/|| 2 < c ie - f (||/|| 2 + ||V/|| 2 ) 



for all t > and some constant C\ . 

\f~L)gds = 



at Jo ' 
L n )) with 



Next, since J?J* cos(.s 
C 1 fR+,L 2 



cos(tVT/)g , we see that j^cos{s\/L)g ds is in 



/ cos(sVL)gds\\ 2 +\\§- t [ cos(sVZ)g ds\\ 2 < c(l + t)e M \\g\ 
Jo Jo 



2 ■ 



Finally, letting v :- 



9 





we obtain J* cos(s\/L)g ds G C°(R + ,V) with 



cos s >/L)</ds|| 2 < \\B 

't 



-ii 



/ cosf sBD)vds\ 



<2\-\\\BD / cos(sBD)vds\\ 2 + \\C\\\\ / cos(sSD)t; tZs|| 2 



< 2A" 1 (|| sm(tBD)v\ 



C\\dte ut \\v\\ 2 ) < c 2 (l + t)e^\\g 



for all t > and some constant c 2 . On using these bounds in f)5.ip . we obtain the 
required estimate. (When B is invertible, then the term (1 + 1) is not needed.) 

Remark 5.9. When a = 1, then L is the operator in L 2 (Q) associated with the 
sesquilinear form J a '■ V x V — > C defined by 

„ n n n 

J A (f,g): = / ^^(^/)^ + Z)Ao/^ + X;A fc (9 fc /)^ + Ao^ 

^ n 7,fc=l 7=1 fe=l 



I 

V 



7 
V 



g) 



for all /, g G V . See, for instance, Chapter IV in 

Remark 5.10. The choice of V = W 1,2 (Q) gives L with Dirichlet boundary con- 
ditions, whilst V = W^ 2 (VL) gives L with Neumann boundary conditions, though 
usually it is assumed in this case that the boundary of Q is at least Lipschitz. When 
V consists of all functions in W 1,2 (Q) which are zero on part of a Lipschitz boundary, 
then the corresponding operator L satisfies mixed boundary value conditions. In [B], 
the authors obtained Davies-Gaffney estimates for the resolvents of such operators 
under similar hypotheses as here, and in a related fashion. 
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